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Detailed high resolution observations of the innermost central parts of nearby galaxies have 
revealed the presence of supermassive black holes 1 " 4 . These black holes may interact with 
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their host galaxies by means of 'feedback' in the form of energy and material jets; this feed- 
back affects the evolution of the host and gives rise to the observed relations between the 
black hole and the host 5 . Here we report observations of the ultraviolet emissions of mas- 
sive early-type galaxies. We derive an empirical relation for a critical black-hole mass (as a 
function of velocity dispersion) above which the outflows from these black holes suppress star 
formation in their hosts by heating and expelling all available cold gas. Supermassive black 
holes are negligible in mass compared to their hosts but nevertheless seem to play a critical 
role in the star formation history of galaxies. 

The near-UV (NUV) detector of the Galaxy Evolution Explorer satellite GALEX 6 covers a 
range in wavelength between 1771 and 283lA and is extremely sensitive to young stellar popu- 
lations. With it, we can detect small mass fractions of 1-3% of young stars formed within the 
last billion years 7 8 . This high sensitivity allows us to trace ongoing residual star formation in 
present day early-type galaxies (ETGs). We select a volume-limited sample of visually inspected 
ETGs from the Sloan Digital Sky Survey 9 that have been observed by GALEX. We remove active 
galactic nucleus (AGN) candidates from our sample to avoid confusion between the observed UV 
emissions from the AGN and from young stars x . 

We use NUV — r colour as a probe of small amounts of recent star formation (RSF) and find 
that the fraction of ETGs that show signs of RSF is strongly correlated with the stellar velocity 
dispersion a, (hereafter "RSF-cr relation"), but not with the luminosity of the galaxy (see Figure 
1). The correlation with a cannot be explained by a systematic trend of a with internal extinction, 
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because such a trend would be exact opposite 11 to that observed. 

We use a semi-analytic model of galaxy evolution in the ACDM paradigm to interpret these 
results 2 - 314 (for model details see Supplements). Our simulation includes all the systematic pro- 
cesses that affect galaxy evolution and takes into account the different merger and star formation 
histories of galaxies 15 . When we run this simulation and investigate the predicted amount of recent 
star formation in low redshift ETGs, we find that all galaxies underwent enough star formation 
within the last billion years to be detected by the extremely sensitive GALEX NUV detector. Mas- 
sive galaxies continued to gain cold gas as fuel for star formation from mergers with satellites and 
through the cooling of hot gas present in those galaxies. 

This overproduction of stars in massive galaxies is a known problem in models of galaxy 
evolution. In order to prevent star formation, a shutdown mechanism is required: cold gas must 
be heated up and expelled to prevent it from turning into stars. Feedback from supernovae (SNe) 
was once a prime suspect, but it turned out to be too low to be effective in massive galaxies 
which have deep potential wells 16 . All the systematics implicit in the ACDM paradigm, including 
feedback from SNe have already been included in our simulation, so we turn to feedback from 
active galactic nuclei (AGN) powered by supermassive black holes (SMBH) as the remaining 
most likely mechanism sufficiently powerful to shut down star formation in massive galaxies to 
reproduce the results of our observations 1719 . 

Energetic outflows caused by gas falling onto the supermassive black holes thought to reside 
at the centres of galaxies can self-regulate the growth of the black hole and the star formation in 
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the galaxy by heating the available gas and blowing it out 5 . Numerical simulations have shown 
that the energy input from a black hole in such an AGN phase is capable of doing this 20 . As SNe 
feedback is incapable of producing the quiescent galaxies we observe with GALEX, we conclude 
that AGN feedback is the most likely process by which the available gas is heated and expelled in 
massive ETGs. 

In addition to this, the nature of the observed correlation with sigma (see Figure 1) allows 
us to infer more about the way in which AGN feedback regulates galaxy evolution. The masses of 
black holes correlate strongly with various properties of the host galaxies and the relation thought 
to have the smallest intrinsic scatter is that between M, and the velocity dispersion a as a proxy 
for the depth of the gravitational potential 1,3 ' 4 . By considering the maximum strength of an active 
galactic nucleus given in terms of the Eddington luminosity L edd oc M, and expressing the gravita- 
tional potential depth by the stellar velocity dispersion a, it is possible to derive a minimum black 
hole mass for which the outflow is capable of depleting all gas from the potential 5 . According to 
our observations, massive ETGs which have large values of a, are on average less likely to have 
undergone recent star formation compared to less massive ETGs: the "RSF-cr relation". Thus, for 
a given velocity dispersion, we assume that there is a critical supermassive black hole mass above 
which the feedback powered by the black hole is sufficiently powerful to heat up and/or expel all 
gas in the host galaxy, thus preventing it from forming any stars even at the 1% level. 

The observed intrinsic scatter in the M,—a relation gives a range of black hole masses for a 
given value of a. Hence, in order to reproduce the decreasing fraction of RSF ETGs with increasing 
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a, the M. jC — a relation should be less steep than the M.—a relation itself. We can derive an 
empirical M #)C — a relation from our observed trend assuming a slope and scatter for the M,—a 
relation. The M, — a relation is usually expressed as a power law of the form M, = a (<r/200) ai , 
so we also try to find such parameters a and ai in the new M #;C — a relation that result in the 
observed "RSF-cr relation". We generate black hole masses for the galaxies in our sample by 
putting their velocity dispersions into the M, — a relation given by Tremaine et al. 3 and perturbing 
the black hole masses by the scatter of the relation; we do this ten times for each galaxy (Figure 2 
(a)). Then we vary a and a\ until we find the best match to the observed RSF-cr relation. We find 
that with 95% confidence M. jC = (1.05±g;^) x lO 8 ^) 3 ' 2 -^. As we are studying a volume and 
magnitude-limited sample, the M. jC — a relation derived from the Monte-Carlo simulation is only 
an approximation and we need to perform full simulations, as they are complete in a as well as in 
absolute magnitude, to derive the actual M. )C — a relation. 

We perform a more realisitc simulation using our semi-analytic models of galaxy evolution 
in order to test the plausibility of the hypothesis of a critical black hole mass controlling star 
formation via AGN feedback. This is a standard semi-analytic simulation in the ACDM paradigm 
that already includes various heating and cooling processes (e.g., SNe). Our model does indeed 
reproduce the observed trend, but only when AGN feedback is included, and the best agreement 
with the observations (Figure 1) is achieved by the relation M. c = 1.26 x 10 8 (^q) 3 - 65 (Figure 2b), 
which is in good agreement with our initial guess. The agreement in the lowest a— bin is poorer 
than in the others due to an increased scatter in our simulated Faber- Jackson relation relative to the 
observed one 21 . The supermassive black holes grown in our simulation have an intrinsic scatter of 
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a factor of 1.5 in mass and are in good agreement with the relation suggested by Tremaine et al. 3 

As our M. jC — a relation is less steep than the M, — a relation we expect an increasing 
fraction of star forming early-type galaxies with decreasing velocity dispersion until eventually 
no AGN will be sufficiently powerful to shut down all star formation. Boselli et al. 4 find that the 
NUV-V (optical) and NUV-H (near IR) colour-magnitude relations for the early-type population of 
the Virgo cluster shows precisely this kind of trend: low luminosity - and therefore less massive - 
Virgo early-type galaxies are increasingly bluer and eventually all lie below the RSF cut-off. We 
use stellar models 23 to convert our RSF criterion to the NUV-V colours of the Virgo galaxies and 
find that below a ~ 80kms _1 all early-types are classified as RSF. Our simulations agree with this 
cutoff very well. 

The possible implications of the M. jC — a relation on galaxy formation and the growth of 
black holes can be summarized by three main regions (see Figure 3): First, the low c-region, in 
which the majority of black holes are too small to be able to halt star formation in their host galaxies 
by feedback during their initial accretion phase. During this phase, black holes grow mainly by 
accreting gas at the Eddington rate and are only limited in their growth by the amount of available 
fuel. Star formation on the other hand is regulated by the energy output from supernovae into the 
interstellar medium and will be proportional to the surface density of the gas 24 . This region extends 
up to a ~ 80km s _1 as suggested earlier. The next phase is at 80kms _1 < a < 240kms~ 1 . Black 
holes are now becoming on average massive enough to be able to terminate star formation by their 
feedback and at the same time limit their available fuel. The mass growth of the black holes is still 
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dominated by accretion but mergers between black holes start playing an increasingly important 
role. Most host galaxies will be identified as AGN or QSO when observed during the accretion 
phase onto the black hole. We sub-divide this region into two: at the lower end, a minority of 
galaxies will have critical black holes; the black hole still grow mostly via accretion. At the high 
end, most black holes are critical and so most of the mass growth now occurs via mergers. The 
boundary between these two is at the intersection between the M. — a and the critical relation, 
i.e. a ~ 165kms _1 . Finally, for galaxies with a > 240kms~ 1 the vast majority of black holes 
is massive enough to halt all star formation and so star formation becomes virtually prohibited 
in such massive galaxies. Both the galaxy and the supermassive black holes now mainly grow 
via dry mergers 25 . Further observations of nearby galaxies hosting supermassive black holes in 
combination with star formation estimates will help to shed more light into the scenario we propose 
and will place a strong test for the existence of "a critical black hole mass" in galaxies. 
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Figure 1 : The relationship between stellar velocity dispersion, luminosity and the fraction of 
star-forming early-type galaxies In panels (a) and (b), we show the observed UV colour-sigma 
and colour-magnitude relations. In the case of the colour-sigma relation, the fraction of blue early- 
type galaxies declines sharply as a function of sigma. In panels (c) and (d), we show the derived 
fractions of RSF (recent star formation) galaxies as a function of velocity dispersion and absolute 
magnitude. In red, we show the observed fractions, which include those early-type galaxies so 
faint in the UV that they are not detected by the GALEX NUV detector. The blue points are the 
fractions of those galaxies in the model that correspond to having had between 1 and 3% of stellar 
mass growth within the last 1 billion years. Error bars show Poisson errors. The vertical lines 
delineate the equal-number a bins. 
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Figure 2: The M,—a and the critical supermassive black hole mass-cr relation (a) The observed 
M, — a data points (black dots) with their errors taken from Tremaine et al. 3 The solid line is 
the best fit to the M. — a relation, with the dash-dotted line showing the l-a scatter. Based 
on this relation we perform a Monte-Carlo simulation for our galaxies (small dots). We colour 
these simulated galaxies depending on whether they lie above or below the critical relation in 
red (quiescent) and blue (RSF) that we derived by fitting the RSF fraction to the observed values 
(Figure 1). The vertical lines delineate the equal-number a bins of the data, (b) We show our 
semi-analytic galaxy models adopting M. c — a relation shown in Figure 2(b), in comparison to 
the empirical M, — a relation. Red and blue dots are the model galaxies that did and did not have 
recent star formation in the last billion years. A new demarcation slope that effectively separates 
these galaxies has a a 3 65 dependence (dashed green) which is close to the Monte-Carlo one in (a). 
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Figure 3: A schematic view of how the critical supermassive black hole mass - velocity disper- 
sion relation regulates galaxy evolution The black hole masses and velocity dispersions of our 
simulated galaxies down to a ~ 50kms _1 . The full points are early-type galaxies, while empty cir- 
cles are late-type galaxies. Red and blue symbols denote quiescent and RSF galaxies, respectively. 
We show the observed M. — o relation (solid line) together with the best SAM critical relation 
(dashed). We indicate the various regimes by colour. 
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Supplements 



1 Galaxy Formation Modelling 

The main strategy behind the modelling approach we follow is to first calculate the collapse and 
merging history of individual dark matter halos, which is governed purely by gravitational interac- 
tions, and secondly to calculate the more complex physics of the baryons inside these dark matter 
halos, including e.g. radiative cooling of the gas, star formation and feedback from supernovae by 
simplified prescriptions on top of the dark matter evolution. Each of the dark matter halos consist 
of three main components which are distributed among individual galaxies inside them. A stel- 
lar, cold and hot gas component, where the latter is only attributed to central galaxies, the most 
massive galaxies inside individual halos. In the following sections, we will describe briefly the 
recipes used to calculate these different components which are mainly based on recipes presented 
in Kauffmann et al. 5 (hereafter, K99) and Springel et al. 3 (hereafter, SOI), and we refer readers for 
more details on model implementations to their work and references therein. 

Throughout this paper we use the following set of cosmological parameters: f2 = 0.3, 
tt A = 0.7, Q b /tt = 0.15, a 8 = 0.9 and h = 0.65. 

Dark Matter Evolution We calculate the merging history of dark matter halos according to the 
prescription presented in Somerville et al. 6 . This approach has been shown to produce merging 
histories and progenitor distributions in reasonable agreement with result from N-body simulations 
of cold dark matter structure formation in a cosmological context 7 . The merging history of dark 
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matter halos is reconstructed by breaking each halo up into progenitors above a limiting minimum 
progenitor mass M min . This mass cut needs to be chosen carefully as it ensures that the right galaxy 
population and merging histories are produced within the model. Progenitor halos with masses 
below M min are declared as accretion events and their histories are not followed further back in 
time. Progenitors labelled as accretion events should ideally not host any significant galaxies in 
them and be composed mainly of primordial hot gas at the progenitor halo's virial temperature. To 
achieve a good compromise between accuracy and computational time, we use M min = 10 10 M Q , 
which ensures that the results presented here are unaffected by numerical resolution effects. 

Baryonic Physics As mentioned above, once the merging history of the dark matter component 
has been calculated, it is possible to follow the evolution of the baryonic content in these halos 
forward in time. We assume each halo consists of three components: hot gas, cold gas and stars, 
where the latter two components can be distributed among individual galaxies inside a single dark 
matter halo. The stellar components of each galaxy are further divided into bulge and disc, to allow 
morphological classifications of model galaxies. In the following, we describe how the evolution 
of each component is calculated. 

Gas Cooling & Reionisation Each branch of the merger tree starts at a progenitor mass of M min 
and ends at a redshift of z = 0. Initially, each halo is occupied by hot primordial gas which 
was captured in the potential well of the halo and shock-heated to its virial temperature T vir = 
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V c / (km s x ) K, where V c is the circular velocity of the halo 2 5 . Subsequently this hot gas 



component is allowed to radiatively cool and settles down into a rotationally supported gas disc at 



15 



the centre of the halo, which we identify as the central galaxy 2 8 - 9 . The rate at which hot gas cools 
down is estimated by calculating the cooling radius inside the halo using the cooling functions 
provided by 10 and the prescription in SOI. In the case of a merger between halos we assume 
that all of the hot gas present in the progenitors gets shock-heated to the virial temperature of the 
remnant halo, and this gas can cool down only onto the new central galaxy which is the central 
galaxy of the most massive progenitor halo. The central galaxy of the less massive halo will 
become a satellite galaxy orbiting inside the remnant halo. In this way, a halo can host multiple 
satellite galaxies, depending on the merging history of the halo, but will always only host one 
central galaxy onto which gas can cool. The cold gas content in satellite galaxies is given by the 
amount present when they first became satellite galaxies and does not increase, but decreases due 
to ongoing star formation and supernova feedback. 

In the simplified picture adopted above, the amount of gas available to cool down is lim- 
ited only by the universal baryon fraction Vl b h 2 = 0.0224 11 . However, in the presence of a pho- 
toionising background the fraction of baryons captured in halos is reduced 12, 13 and we use the 
recipe of Somerville et al. 14 , which is based on a fitting formulae derived from hydrodynamical 
simulations 13 , to estimate the amount of baryons in each halo. For the epoch of reionisation, we 
assume z reion = 20, which is in agreement with observations of the temperature-polarisation cor- 
relation of the cosmic microwave background 15 . 

Star Formation in Discs and Supernova Feedback Once cooled gas has settled down in a disc, 
we allow for fragmentation and subsequent star formation according to a parameterised global 



16 



Schmidt-Kennicutt law 16 of the form M* = aM co i d /t dynjgal , where a is a free parameter describing 
the efficiency of the conversion of cold gas into stars, and t dyn>ga i is assumed to be the dynamical 
time of the galaxy and is approximated to be 0.1 times the dynamical time of the dark matter halo 5 . 
As in K99 we allow star formation only in halos of V c < 350 km/s to avoid excessively-massive 
central galaxies in clusters. 

Feedback from supernovae plays an important role in regulating star formation in small mass 
halos and in preventing excessively-massive satellite galaxies from forming. We implement feed- 
back based on the prescription presented in K99 with 

v , , 4 rj SN E SN 

AM reheat = -e — — — AM*. (1) 

3 V c 2 

Here we introduce a second free parameter e which represents our ignorance about the efficiency 
with which the energy from supernovae reheats the cold gas. The expected number of supernovae 
per solar mass of stars formed is given by tjsn = 5 x 1CT 3 , taken as the value for the Scalo initial 
mass function 17 , and E S n = 10 51 erg is the energy output from each supernova. We take V c as the 
circular velocity of the halo in which the galaxy was last present as a central galaxy. 

Galaxy Mergers We allow for mergers between galaxies residing in a single halo. As mentioned 
earlier, each halo is occupied by one central galaxy and a number of satellite galaxies depending 
on the past merging history of the halo. Whenever two halos merge, the galaxies inside them are 
going to merge on a time-scale which we calculate by estimating the time it would take the satellite 
to reach the centre of the halo under the effects of dynamical friction. Satellites are assumed to 
merge only with central galaxies and we set up their orbits in the halo according to the prescription 
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of K99, modified to use the Coulomb logarithm approximation of SOI. 

If the mass ratio between the two merging galaxies is M ga i :1 /M galj2 < 3.5 (M galjl > M ga ^ 2 ) 
we declare the event as a major merger and the remnant will be an elliptical galaxy and the stellar 
components and the gas will be treated according to the prescriptions below. In the case of minor 
merger M ga i^/M ga i j2 > 3.5 the cold gas in the disc of the smaller progenitor is assumed to settle 
down in the gas disc of the remnant and its stars contribute to the bulge component of the remnant 5 . 

Formation of Ellipticals and Bulges It is suggested that major mergers will lead to the formation 
of elliptical galaxies 18 . Indeed detailed numerical simulations in the last decade seem to support 
this hypothesis 19 ' 20 , and we wil assume in the following that major mergers disrupt the discs in 
progenitor galaxies as seen in various numerical simulations and relax to a spheroidal distribution. 
During the merger, any cold gas in the discs of the progenitor galaxies is assumed to be funnelled 
into the centre of the remnant where it ignites a starburst which transforms all of the cold gas into 
stars contributing to the spheroidal component 3 ' 5 . The second assumption is certainly a simplifica- 
tion of what might happen since we neglect the possibility that not all of the cold gas is funnelled 
to the centre but some fraction of it may for example settle down onto an extended disk which 
continues to grow inside-out by fresh supply of gas from tidal tails 2124 . The results of Barnes 24 in- 
dicate that 40% — 80% of the initial gas in the disc could end up in the central region of the remnant 
and be consumed in a starburst. The exact amount is somewhat dependent on the merger geometry 
and on the mass ratio of the merger. Unfortunately, a sufficiently large survey investigating the gas 
inflow to the centres of merger remnants is not available yet, so we use the simplified approach of 
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assuming that all cold gas gets used up in the central starburst. This prescription for the fate of the 
cold gas results in an overestimate of the spheroid masses and an underestimate of the secondary 
disc components in our model. This is not very significant for massive ellipticals since they are 
mainly formed in relatively non-dissipative mergers 25 . 

Formation and Growth of Super-Massive Black Holes We here follow the model introduced 
in Kauffmann & Haehnelt 26 in which super-massive black holes get formed and fed during major 
mergers of glaxies. The assumption is to say that a fraction of the the cold gas available in the 
progenitor discs, that is funnelled into the centre of the remnant, will be accreted onto the black 
hole. Kaufmann & Haehnelt introduced the following scaling law for the effectivness of this 
procces: 

»V _ fbhM co id 

*~ 1 + (280km s-VK) 2 ' 
where M co i d is the amount of cold gas available in the progenitor disks, V c the circular velocity 

of the dark matter halo and f bh a free parameter. This scaling provides a good fit to the observed 

relation between the mass of the super-massive black hole and the velocity dispersion of a galaxy 26 . 

Following Croton et al. 18 we here assume the velocity dispersion of a galaxy to be identical to the 

circular velocity of the dark halo the galaxy is embedded in. This however is not strictly true as e.g. 

in the case of an isothermal sphere V c /a = y/2. Several authors have investigated the correlation 

between V c and a- 27 - 28 and found different correlations. We here note that using a different scaling 

between V c and a in our models still produces the same M. — a— relation if we adjust the free 

parameter f bh accordingly. 
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Feedback from Super-Massive-Black Holes We here introduce a new prescription for the mod- 
elling of feedback from super-massive black holes, which is based on the results presented in this 
paper. According to our analyses a critical black hole mass exists for a given galaxy velocity dis- 
persion a, at which feedback is so strong that it reheats all the available cold gas and hence prevents 
further star formation. This critical black hole mass is derived from our empirical fit to the residual 
star formation fraction found in our data as 



In galaxies with black holes more massive than the corresponding critical black hole mass we stop 
any cooling of gas and star formation. We callibrate this relation to give the right fraction of RSF 
galaxies as seen in Figure 1 . 

Free Parameters We normalise our two model parameters for the star formation efficiency a 
and supernova feedback efficiency rj by matching the /—band Tully-Fisher relation of Giovanelli 
et al. 29 and requiring that spiral central galaxies of halos with circular velocity Vc = 220 km/s 
have on average 10 11 M Q of stars and 10 9 M of cold gas 30 . The third free parameter f bh , which 
regulates the black hole growth is set to match the observed M. — ex— relation and is f bh = 0.02. 

2 The Virgo Cluster Data 

We derive the NUV — V colour-c relation for the Virgo cluster down to much smaller values of a 
than in our GALEX-SDSS sample by combining the near-UV and optical photometry from Boselli 
et al. 4 with velocity dispersion measurements from the GOLDmine database 31 . The purpose of 




(3) 
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looking at the Virgo cluster galaxies is to constrain the low end of the a and so get an idea of the 
point at which all galaxies are RSF galaxies and thus below critical. Since the optical photometry 
of the Virgo galaxies in this sample are in V-band, we approximately convert the RSF criterion of 
NUV-r=5A to NUV-V. We use simple stellar populations of solar metallicity and three ages (3,6 
and 9 Gyr) to derive the NUV-V RSF criterion; these are NUV-V= 5.16, 5.12 and 5.09 respectively. 
In all three cases, the point at which no Virgo early-type galaxy is above this cutoff - i.e. all show 
some signs of young populations - is at a ~ SOkms^ 1 . 



3 The RSF criterion 



We use NUV-r colour as a discriminant between those galaxies which are quiescent and those 
which show signs of recent star formation. Besides AGN, there is one further effect which can 
mimic young stars in the near-UV. Extremely old stellar populations can give rise to the UV upturn 
phenomenon in some early-type galaxies 32 33 . We choose the NUV-r colour of one of the strongest 
UV upturn galaxies NGC 4552 33 to be the boundary between galaxies with no young component 
on the one hand and those that are so blue they must have some young stars in them. Although 
still quite limited, theoretical population synthesis models which are supported by empirical data 
(Lee et al. 34 ) also suggest that the NUV-r of passively evolving ETGs could not be much bluer than 
5.4. Using this criterion, we compute the fraction of galaxies that must have experienced some 
star formation within the last billion years. The ages and mass fractions of these young stellar 
components are generally 300-500 Myr and 1-3% by mass 35 . Due to the high sensitivity to dust 
extinction, these fractions are a lower limit as cold gas and dust have been detected in many early- 
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type galaxies 36 37 . NUV-r probes recent star formation for up to 1 Gyr and so traces different time 
scales than those of AGN feedback, so there may be a certain amount of time lag between the UV 
emissions and the processes underlying feedback. 

4 Visual Inspection of Galaxy Morphology 

The process and criteria of the visual inspection are described in detail in Schawinski et al. 1 . In 
summary, we find early-type galaxy (ETG) candidates from the SDSS database using fracJJev 
parameter, which is the weight of the de Vaucouleurs component in the two-component (de Vau- 
couleurs and exponential disc) fits. We use fracJJev > 0.95, which is a highly conservative 
criterion, hoping to exclude as many spiral interpolers as possible. Despite this, still some spi- 
ral bulges remain in our sample. Hence, we remove them via visual inspection on gri bands 38 . 
Since visual inspections are also subject to errors, we investigat out to what redshift and apparent 
magnitude visual inspection based on SDSS images is reliable by comparing them with one of 
the COMBO- 17 fields overlapping SDSS. The COMBO- 17 image is significantly deeper and has 
much better seeing with 0.7" as opposed to 1.4" which are typical for SDSS. By comparing our 
classification based on these two different data sets, we concluded that visual inspection of mor- 
phology using SDSS was possible to z ~ 0.13 and R < 17.31. We settle on more conservative 
limits of z < 0.1 and r < 16.8. The presence of late-type interlopers with star-forming disks or 
spiral arms is therefore not a significant concern and does not give rise to the correlations between 
M r , a and NUV-r color. A selection of sample images with their velocity dispersions are shown in 
Figure |4j 
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Figure 4: Example galaxies from our SDSS sample. The images are gri composites 3 and span a 
wide range in velocity dispersion. 
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